Multiple-scale approach for the expansion scaling of superfluid quantum gases 
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We present a general method, based on a multiple-scale approach, for deriving the perturbative 
solutions of the scaling equations governing the expansion of superfluid ultracold quantum gases 
released from elongated harmonic traps. We discuss how to treat the secular terms appearing 
in the usual naive expansion in the trap asymmetry parameter e, and calculate the next-to-leading 
correction for the asymptotic aspect ratio, with significant improvement over the previous proposals. 
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I. INTRODUCTION 

The physics of ultracold quantum gases has been the 
object of an intense experimental and theoretical inves- 
tigation since the achievement of Bose-Einstein conden- 
sation in dilute trapped gases in 1995 [H [2] . A common 
ingredient of these experiments is the so called time-of- 
flight, that is the expansion of the gas after the release 
of the confining potential. For some typical experimental 
regimes, the expansion after the release from harmonic 
traps can be described by scaling solutions, in terms of 
three scaling parameters bi{t) {i = 1,2,3), which obey a 
system of second-order ordinary differential equations, of 
the general form (where uji are the trapping frequencies) 
[IH6] 
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with initial conditions 6^(0) = 1, 6i(0) = . In particular, 
7 = 1 corresponds to the case of a Bose-Einstein conden- 
sate (BEC) in the Thomas-Fermi (TF) hmit [2 H i], 
whereas 7 = 2/3 refers to a unitary superfluid Fermi gas 
[2 O E] . The solution of the above equations allows the 
complete characterization of the expansion of the cloud 
size as Ri{t) — Ri{0)bi{t). 

In the case of strongly elongated cylindrically- 
symmetric traps, e — uJz/uia_ <C 1, the above equations 
read 
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(3) 



with T = uj±^t. Usually, these equations are solved by 
means of an expansion in powers of e^ , by the ansatz 
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b,{r) = Y,e^-t\r)- 
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to the well-known Castin and Duni scaling [3] 

b^{T) = v/l + r2 (5) 

bz{T) = l + (^ (t arctan t - In \/l + tA (6) 

that nicely describes the experimental data [T]. 

From these expression one can also infer the asymp- 
totic aspect ratio, defined as \vLnT^^{Rj_(T) j Rz(t)) — 
(-R_L(0)/i?z(0))(2/7re2) [T]. However, one may notice that 
the above expansion is strictly valid only in the pertur- 
bative regime, when the second order correction in ([6]) is 
smaller that unity, namely for r ^ 2/'r:t^. This condi- 
tion is in general well satisfied in the typical experimental 
regimes, but in principle does not permit the direct ex- 
traction of the asymptotic limit, because of the secular 
term r arctan r that eventually invalidates the perturba- 
tive expansion. 

Here we present a general method, based on a multiple- 
scale perturbative expansion, that allows the derivation 
of a uniformly valid expansion, where the hierarchy of 
sub-leading terms is preserved at any time. This can be 
achieved by means of a proper resummation of the secular 
terms [7]. 

The paper is organized as follows. In Sect. [IT] we re- 
formulate the problem by using a Hamiltonian approach 
(see also [4]). In QIA we show that the leading term of 
the asymptotic ratio can be computed already at zeroth- 
order in the naive expansion, by using the expression in 
terms of the canonical momenta (avoiding secularities) . 
We find that this result may represent a good approxi- 
mation of the exact (numerical) asymptotic aspect ratio, 
depending on the value of 7, and provided that e is small 
enough. However, it may deviate significantly from the 
exact solution for 7 < 1, alread y for not too large val- 
ues of e (e.g. e « 0.05). In QIB we show that the 



For example, for the BEC case 7 = 1, one retains the 
zeroth-order in &j^, and up to second order in 62, leading 



next-to-leading order in the naive expansion gives rise to 
secular terms for any 7, requiring therefore a different 
expansion approach. Then, in Sect. |III| we introduce 
the general ideas for the multiple scale approach, and 
consider the next-to-leading correction to the asymptotic 
aspect ratio, addressing in particular the case 7 = 1 and 
1/2 < 7 < 1. In both cases the expansion tuns out to be 
non-analytic due to the presence of terms proportional to 



e^ loge (7 = 1) or e^^^''' ^) (7 < 1). We then recapitulate 
and offer prospective applications of the method. 



II. FORMULATION OF THE PROBLEM 

Equations (|2|-(l3| are equivalent to the canonical equa- 
tions {bi ~ dH/dpi, Pi = —dH/dhi) 
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generated by the following Hamiltonian |3] 
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where the initial conditions &i(0) = 1, 6^(0) — set 
the total energy to H = I/27. Asymptotically, 6^ « 
and therefore bi ex r, so that Ai?oo = limr->.oo(&j_/^z) = 



The asymptotic value of pk is given by the integral 
lim p(°) = / dTp(0)(^) 



/7dx 
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with 7 > 1/2; here we have used the fact that pi = 
bj_ pj_ /2, combined with Eq. ( 13 ). From the same equa- 
tion we get 



lim p^j^' 
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Therefore, within this hamiltonian formalism, the leading 
term of the asymptotic ratio can be computed already 
at zeroth-order by using the expression in terms of the 
canonical momenta, ARao — {^/2€'^)ipj_/Pz)cx> 9J. This 
has to be compared with the usual approach in terms of 
the bi, where one should include terms up to order e^ 
in bz (see Eq. ([s]) and [3]). By using the above leading 
values of the momenta we obtain 
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A. Zeroth-order of the naive expansion 



From this expression we can extract the TF bosonic case 
7=1, and the unitary fermionic case, 7 = 2/3 



Let us start by considering the naive expansion in Eq. 
Q, both for bi{T) and Piir). To lowest order in e^, the 
Hamilton equations are 



^f-pf pf 



6i°) = p(°) = 



1 



^(0)^.+i^(0). 



2;,(o)^-'^(o)-'+^ 



(10) 
(11) 



Then, the zeroth-order solution is simply bi — 1, and 
the equation for b^ becomes 
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From energy conservation we have 
1 
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and b: may be easily obtained in implicit form as 
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Here we have defined for latter use 
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that coincide with the results in (TJ [S] . 

In Fig. [Tlwe compare this zeroth-order prediction (dot- 
dashed line) with the asymptotic aspect ratio obtained 
from the full numerical solution [10 of Eqs. (|2])-(l3]) (solid 
line), as a function of 7. This figure shows that the 
zeroth-order result (17 1 is rather good for 7 > 1, but de- 



viates significantly from the exact (numerical) solution 
for lower values of 7, when e is not sufficiently small. We 
will see in the following how to improve this result. 

Notice also that we have not yet implemented the can- 
cellation of secularities (that in fact occurs at the next 
order in the expansion). 



B. Next-to-leading order (naive) 



The next-to-leading order equations are 
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FIG. 1. (Color online) Asymptotic aspect ratio 2e^ARoo as a function of 7, for e — 0.01,0.05,0.1 (from left to right). The 
full numerical solution of Eqs. ([2|-([3| (orange solid line) is compared with the zeroth-order expression in Eq. ( |17[ ) (green 
dot-dashed line) and with that including the next-to-leading corrections (magenta dashed line, red dot at 7 = 1) discussed in 
Sect. Ill . The vertical (dashed) lines correspond to 7 = 2/3 and 7=1. 



Let us consider Eq. (22 1; by using Eq. (15) the asymp- 

totics is 6z ^ ^(7), and the corresponding large time 
behavior of the axial scaling parameter 



implicit expression 



b,{T)r^l + e^c{-f)T, 
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grows unboundedly with r. This behavior is physically 
correct because we expect non-zero asymptotic limits for 
both velocities. However, with regard to perturbation 
theory, the growing of the first order correction must be 
interpreted as the occurrence of a secular term relative 
to the naive zeroth-order solution bi = 1. That is, the 
perturbative hypothesis breaks down after times of order 
l/e^c(7). Therefore, to avoid the appearance of these 
terms we propose a reinterpretation of the perturbative 
analysis in the framework of multiple-scale analysis. It is 
notable that the resummation involved in this procedure 
implies that the higher-order terms of the velocities are 
O ((e^)^''"^) if 1/2 < 7 < 1, and O(e'^lne) when 7 = 1. 



III. MULTIPLE SCALE APPROACH 

We begin by introducing an additional time scale 
Ti(t) = e^r, and assume a perturbative expansion 

6,(r,ri) = Qf )(r,ri) + e^Q\'\T,n) + ... (25) 



K(T,ri)=/^<"'(r,Ti) + e2/^(i)(r,ri)-)- 



(26) 



in such a form that the corrections Q] and P^ must be 

not secular with respect to Qi and i^ in r time; i.e. 

Q\ and P^ must be consistently be of order e^ with 

respect to Ql and Pj for all values of r. The derivative 
symbol "•" is now replaced by d/dr = dr + e^dn- At 
leading order, it is nearly obvious that the choosing of 

)i°)(r,n) = gf(ri)^l + c(7)n, (27) 
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removes the linear secularity. An explicit computation 
with the assumed initial conditions yields the following 
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which produces the asymptotic behavior 



Pf(r,rr) 



Pi"Hr,r,) 



V7 




7^-^(«i°^(n))'%^-^-^ 


2-47 


1 

V7 


i-^7^-^(9f(n))'V-^- 


c(7) 


^7-(gi°^(n))"^%-^- 



(28) 



(29) 



(30) 



(31) 



2 2-47 

for r — > 00 at fixed ti (ti = 0(1)). Now the equation for 
Qz is simply 



dQ 



(1) 



dr 



-c(7) + 2Pi°) 



(32) 



which in the region r — > 00, ti = 0(1) is of order t^^^'*'. 
This leads to an asymptotic growth rate that is slower 
than linear. 

Let us now consider the cases 7=1 and 1/2 < 7 < 1 
separately. For the boson case 7 = 1 the expressions are 
more transparent and become 



Qf(T,Ti)- 



^f(r,n) = 
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Therefore, we see that Qz (r, ti) only grows logarithmi- 
cally as r — > oo and ri is 0(1), in contrast with the naive 
perturbative result of Eq. ([6| . 

We can exploit these results to improve the perturba- 
tive evaluation of the aspect ratio to higher order in e. 
Notice that earlier we were able to predict the asymp- 
totic ratio to first order from the zeroth order result by 
using the Hamiltonian approach; analogously here we can 
obtain an improved aspect ratio. 

Let us now consider the axial momentum, whose 
asymptotic value can be obtained as 
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The corresponding integrations produce 
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(42) 
and 
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that corresponds to the first order correction to the 
asymptotic radial momentum. 

Thus, one obtains the asymptotic values 



in order to extract the next-to-leading corrections, it is 
sufficient to consider its equation of motion up to order 
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Then, the integration of the first term gives 
dr 
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where the "resummation" has been crucial to ensure the 
convergence of the integral. For the other two terms, 
we may neglect the dependence on ti (since it intro- 
duces corrections of order e^ in the integral, which com- 
pete with higher order terms in the corrected pertur- 
bative expansion), by using qi (0) = 1, Qi (r, 0) — 
rarctauT — In \/l + t^ — |t, and Q\_ {t,0) 



VT 



Therefore Q^ ' (r, 0) = bj_ (r) is exactly the first non- 
rcsummed radial correction satisfying 
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where bi (t) = t arctan t — In vT+^ (see Eq. (|6| . The 



solution of Eq. ( 39 ) with zero initial conditions may be 
written as 



pT = i- 
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(44) 



pT^^ + ^ (2 + 31n2 + 41n(7re2)) +0(e4in£)(45) 

These results represent a significative improvement for 
the asymptotic ratio ARao = (l/2e^)(p_L/pz)oo with re- 
spect to the zeroth-order result in Eq. (18), as shown in 
Fig. [2] (red dotted line) , and in the three panels of Fig. 
[1] (red dots). 

Finally, let us now consider briefly the asymptotics of 
momenta in the case 1/2 < 7 < 1. Now pf^ is given by 
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where the O(e^) terms come from integrations that are 
analogous to those in the last brackets of Eq. (38). It 



is possible to perform an asymptotic expansion of above 
integral and to obtain the two leading terms. The result 
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with G{t, s) being the following Green's function 
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FIG. 2. (Color online) Asymptotic aspect ratio 2e^j4_Roo as 
a function of e^, for 7=1. The full numerical solution of 
Eqs. ([2|-([3| (orange solid line) is compared with the zeroth- 
order expression in Eq. ( 17 1 (dot-dashed line) and with those 
including next-to-leading corrections for 7 = 1 (red dotted 
line). 



pT-^+ 6(1 27) ' [^^ '^'^^ + 7, 1 - 27, 2 - 27,-1) + (1 - 27) 2Fi(l + 7, 37, 1 + 37, -1)] 
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Notice that the expansion contains a term proportional 
to (e^)^'''"^ that dominates over the e^ term when 1/2 < 
7 < 1. However, the latter term cannot be discarded, as 
it is needed to cancel the singularity in the limit 7 — > 1, 
producing an O(e^lne) term. Therefore, the improved 
asymptotic ratio ARao can be obtained by combining Eq. 
07\ with 



FIG. 3. (Color online) Asymptotic aspect ratio 2e^ARcx, as 
a function of e^, for 7 — 2/3. The full numerical solution of 
Eqs. ([2|-(|3| (orange solid line) is compared with the zeroth- 
order expression in Eq. (171 (dot-dashed line) and with that 
including next-to- leading corrections (dashed hne). 
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where the second term is easily obtained from energy 
conservation. Even in this case, it represents a signifi- 
cant correction with respect to the zeroth-order result, 
as shown in Fig. [T] 

In particular, the asymptotic expressions for the 
fermionic case, 7 — 2/3, are 



p^ ^ 1.6057^ 
p"^ = 1.2248- 



2.5333 6^/3 + 
2.1051 e^-h.. 



(49) 
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and the corresponding aspect ratio is shown in in Fig. 
^ as a function of e^. In this case, the improvement 
provided by the multiple-scale approach over the naive 
zeroth-order result is quite significant in the small e 
regime. 
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IV. CONCLUSIONS 

The asymptotic aspect ratio predicted from the naive 
expansion is surprisingly adequate, when considering 
that the naive expansion is invalid for the large times 
for which the aspect ratio is desired. We have explained 
this success as being derived from the Hamiltonian char- 
acter of the equations of motion, so that in fact it is not 
a perturbative result. We have also produced resummed 
corrections to next order, in which the non-analycity be- 
comes apparent. The techniques applied here have a long 
history in applied mathematics; most of the published ex- 
amples of the multiple scales method, however, concern 
bounded periodic motion. Here we make an novel ap- 
plication thereof to asymptotically linear solutions. We 
have provided perturbative expansions in the anisotropy 
parameter which are valid over the whole expansion time, 
and not only for small times. This could be used to ex- 
tract trap information out of longer time-of-flight exper- 
iments. The techniques presented here have wider appli- 
cability, and could be applied to, for instance, expansions 
of multiple species. 
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